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We study the phase diagram of quark matter at finite temperature and density in the 
strong coupling lattice QCD with one species of unrooted staggered fermions including finite 
| coupling (l/g 2 ) effects for color SU(A r c ). We find that we may have partially chiral re- 

stored medium density matter at N c = 3, which would correspond to the quarkyonic matter 
suggested at large N c . 



The QCD phase transition at finite temperature (T) is the latest vacuum phase 
{Sj ■ transition in our universe, and it can be experimentally investigated at RHIC and 

LHC. Phase transitions at finite chemical potential (//) in dense matter may be also 
realized during black hole formation or in neutron stars, where we have the following 
central question. "What is the next to the hadronic Nambu-Goldstone phase in the 
larger (i direction ?" Monte-Carlo (MC) simulations are not yet reliable in the region 
/i/T > 1 because of the notorious sign problem, then it is necessary to invoke some 
approximations such as the large number of colors (A^) 1 ) or the strong coupling limit 
(SCL) 2 )- 9 ) in order to answer this question in QCD. 

Recently, McLerran and Pisarski have shown that the next phase at large N c 
■ should be the so called quarkyonic phase, in which the colors are confined and the 

baryon density is high. 1 - 1 At large N c , gluon contribution to the pressure ~ 0(A 2 ) is 
larger than those from hadrons ~ 0(1) and quarks ~ 0(N C ), then the deconfinement 
transition temperature is independent of the quark chemical potential /i as far 
as it is moderate fj, ~ 0(1). In the confined region T < Tj, quark number density 
is exponentially suppressed if fi is below the quark mass m q ~ 0(1), but it rapidly 
grows at fj, > m q and soon reaches high density ~ 0(N C ). This dense matter has 
a characteristic feature that it is made of quarks, while only b&ryonic excitations 
are allowed because it is confined. If the quarkyonic phase is the next at N c = 3, it 
5_i ■ may be formed at high densities in compact astrophysical phenomena or in heavy-ion 



collisions at 10-100 A GeV. MC results with the density of state method also show the 
transition to high density phase, 10 ) but its nature is not yet known. Very recently, 
quarkyonic matter is found to exist at N c = 3 in effective models of QCD, 12 )' 13 ) 
while the results may depend on the details. 14 ) It is now very important and urgently 
required to discuss the possibility of the quarkyonic matter phase in QCD for N c = 3. 

The above discussion tells us that the quarkyonic transition, the transition from 
the chiral broken Nambu-Goldstone (NG) phase to the quarkyonic matter, is charac- 
terized by the quark number density. If the quarkyonic matter exists, the quarkyonic 
transition occurs at fj, ~ m q where density rapidly grows as p q = 0(1) — ► 0(N C ), 
and the chiral restoration follows at higher chemical potential as shown in Fig. [TJ In 
this paper, we discuss the possibility of the quarkyonic matter in the strong coupling 
lattice QCD (SC-LQCD), 2 )~ 8 ) which is another powerful tool in studying dense mat- 
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Fig. 1. Schematic phase diagram with the quarkyonic matter. 



ter. We take account of the finite coupling effects in the next-to-leading order (NLO, 
1/g 2 ) , 5 )~ 8 ) and introduce an auxiliary field representing the quark number density 
p q as an order parameter in addition to the chiral condensate a. In a previous work 
on the phase diagram with NLO effects, 7 ) the order parameter representing p q was 
not introduced. As shown later, the multi-order parameter (cr and p q ) treatment is 
essential in understanding the quarkyonic transition. In SC-LQCD, we consider the 
situation where g 2 is large while N c is fixed. This condition is somewhat different 
from that assumed in the large N c argument, where N c is assumed to be large and 
the 't Hooft coupling N c g 2 is fixed. As we show later, we find the quarkyonic matter 
to appear in some region of 1/g 2 at -/V c = 3. This observation together with the large 
N c discussion 1 ) support the existence of the quarkyonic matter in a wide region of 
the (N c , l/g 2 ) plane. 

We start from the lattice action with one species (unrooted, four flavors) of 
staggered fermions (x) in the lattice unit, 

S LQC D = ^E[ y+ (' u )-^(' u )]+E m o M - + 5(s) -4E tr [ C7n + C/ n] ' (1) 

X X □ 

where mo and \i are the bare quark mass and the lattice chemical potential, respec- 
tively. The mesonic composites are defined as M x = XxXx-, V x + = e>1 XxUo,xX x+ d 
and V~ = e.~^x x+ (pl x Xx, where sum over color indices are assumed. The spatial 
hopping action of quarks is given as, 

1 d 

S(S) ^EE lxxUj, x x x+ j ~ X x+ jUl x Xx\ , (2) 

X j = l 

where d is the spatial dimension and r)j tX = (— l) Xo_l ^i- 1 is the staggered sign 
factor. The gluon degrees of freedom are described by the temporal and spatial link 
variables and plaquettes, £/o, Uj, U\j. This lattice action is invariant under the chiral 
transformation \ ~^ eldex X with the 75-related factor e x = (— I) Xo_1 ^ Xd . 

In the strong coupling region (g S> I), we can evaluate the plaquette effects 
through the expansion in the power series of 1/g 2 (strong coupling expansion). In 
each order of 1/g 2 or the number of plaquette, we can exactly carry out the integral 
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over link variables by using the SU(A r c ) group integral formulae, f dU U a bUl d = 
$ad3bc/N c and so on, and obtain the effective action of hadronic composites. 2 ) 

Before discussing the NLO effects, we briefly summarize the procedure to ob- 
tain the effective potential in the leading order, i.e. the strong coupling limit (SCL), 
where the sophisticated framework based on the 1/d expansion, n ) mean field approx- 
imation and finite T treatments of the quark deter minant, 3 )' 6 ) has been established. 
After the spatial link (Uj) integral, we obtain the hadronic hopping action shown in 
the third graph in Fig. [2p from the spatial quark action S^ s \ The SCL effective 
action is given as, 

Sscl = \ £ [V + (») - V-M] + £mo^. - | M * M * + j + 



(3) 



where b a = d/2N c . In this work, we adopt the leading order terms of the 1/d 
expansion 11 ^ and omit higher order terms 0(l/y/d). In the 1/d expansion, the lead- 
ing term M x M x+ j is assumed to remain finite at large d, then the quark fields 

scale as x> X x d~ 1 ^ and higher power terms of quarks are found to be suppressed 
as 0(l/y/d) for N c > 3. Through the Hubbard-Stratonovich (HS) transformation, 
the chirally invariant four-fermi term Sscl ^ s converted to the quark mass term 
b a axXi where the finite chiral condensate a = —(M x ) breaks the chiral symme- 
try spontaneously and generates the quark mass dynamically. We can carry out the 
Gauss integral over the quark fields x an d X m the so-called finite temperature treat- 
ment: The determinant of the temporal quark hopping matrix with the anti-periodic 
boundary condition is obtained by utilizing the Matsubara product technique 3 ) or 
the recursion relation. 6 ) It is possible to carry out the temporal link integral of 
this determinant over Uq in the Polyakov 3 ) or temporal gauge. 6 ) Then the effective 
potential is obtained as, 

•^ff CL = y<r 2 +V,K;/z,r) , (4) 



V q (m q ; fi,T) = -Tlog 



Xn +2 cosh 



L T 



(5) 



y ( x smh[(N c + l)E q (m q )/T] 

XNAmq) = s,n H E q (m q )/T] ' (6) 

where we have regarded the inverse of the temporal extension 1 /N T as temperature T, 
and E q {m q ) = arcsinh(m g ) represents the one dimensional quark excitation energy 
coming from the constituent quark mass m q = b a a + tuq. 

We shall now evaluate NLO (1/g 2 ) correction terms coming from the plaquette 
action in Slqcd- We concentrate on the leading order of the 1/d expansion. By 
integrating out spatial link variables, two types of terms shown in the last two graphs 
in Fig. [2] are found to appear from temporal and spatial plaquettes, 5 )' 6 ) 

AS (r) = A.y ,y+y- y+y- ) (7) 

x,j>0 



4 



Letters 



x + 

o 



x + 




V+(x) V-(x) M„M m+j 

Fig. 2. Effective action terms in the strong coupling limit and NLO (l/g 2 ) corrections. Open 
circles, Filled circles, and arrows show x, X> an d U v , respectively. 



Asi s) = , ^ . V M X M ,~-M .*M , t , , (8) 
where (3 T = d/N 2 g 2 and f3 s = d{d - l)/8N*g 2 . 

_ (t) 

The NLO terms contain the product of different types, such as V + V in AS\ . 
In order to treat these terms, we propose a new mean field technique, named Ex- 
tended Hubbard-Stratonovich (EHS) transformation. Let us consider to evaluate a 
quantity e aAB , where (A, B) and a represent arbitrary composite fields and an pos- 
itive constant, respectively. We can represent e aAB in the form of Gaussian integral 
over two auxiliary fields ((p,4>), 



e *AB 



d v d<pe- a { { *~ {A+B)/2)2+{(t> - i{A - B)/2)2 } +aAB 
dtpdcfie-^-iA+B^^A-B^} (9) 

The integral over the new fields (92, 4>) is approximated by the saddle point value, 
ip = {A + B)/2 and (ft = i{A — B)/2. Specifically in the case where both (A) and 
(B) are real, which applies to the later discussion, the stationary value of 4> becomes 
pure imaginary. Thus we replace eft — > iu and require the stationary condition for 
the real value of u, 



e 



aAB ~a{^ 2 -(A+B)^uj 2 +(A-B)u;} 



. (10) 

stationary 



In the case of A = B, Eq. (jlUj) reduces to the standard HS transformation. We find 
that the e aAB is invariant under the scale transformation, A — ► XA and B — ► A _1 -B. 
This invariance is kept in rhs of Eq. ()10p . since the combinations ip — uj = (A) and 
<p + oj = (B) transform in the same way as A and B, respectively. This means that 
the effective potential is independent of the choice of A. 

We apply the EHS transformation to the NLO terms. For the temporal and spa- 
tial plaquette action terms, ASp and ASg S \ we substitute (a, A, B) = {f5 T /Ad, — , 
and ((3 s /d(d - l),M x M x+ ~ p M x+ft M x+ ~ k+ -j) in Eq. flUty, respectively, and obtain, 

(11) 
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V ' x,0<k<j 

(12) 

We can absorb NLO terms, ASp in Eqs. (jlip and (|12p . in the coefficient 
modification of the SCL effective action, S'scl in Eq. ([3|). We assume that auxiliary 
fields take constant and isotropic values, which are independent from the space-time 

(s) 

x and spatial directions j, k. In AS\ , effects of lo s disappear for constant auxiliary 
fields, and ip s modifies the coefficient of MM terms in S'scl as —b a /2d — * —b a /2d, 
where b a = b a + 2(3 s ip s . Combined with the temporal hopping term in SlqcDj 
the coefficients of are found to be, Z^/2, where Z± = 1 + (3 T {4> T ± uj t ). We 
rewrite these coefficients as Z± = Z x exp(±<fyi), then NLO effective action containing 
fermions is written as, 



s$2o =%e i v+ (v - v~m +E^M x - b f d Y: M x M x+] + o(±± 

x,j>0 ' 



(13) 



where jj, = fj, — 6fA = fj, — log \J Z + / Z_ and Z x = \JZ+Z_ represent shifted chemical 
potential and the wave function renormalization factor. 

Now we can repeat the same procedure as the standard SCL prescription. Bosoniza- 
tion of the four-fermi interaction in Snlo leads to the quark mass term (b a cr+mo)M, 
and the quark integral and the temporal link integral are also analogous to the 
SCL case. The effects of wave functional renormalization factor Z x would be a lit- 
tle bit non-trivial: Let us define the temporal hopping matrix as YLxy Xx^xyXy = 
YlxWx ~ Vx~]/2, then the quark determinant in SCL reads Det [K^ +l(b a a+mo)] . 
With NLO effects, this determinant is replaced with, 



Det 



Z x K {r \ji) + liKa + mo)] = Det [Z x ] Det [k (t) (/i) + l(b a a + m )/Z x ] , (14) 



where Det represents the determinant of the temporal and the color matrix. From the 
factor Det[Z x ] in Eq. (|14p . the additional term — N c logZ x appears in the effective 
potential. The second factor in the second line of Eq. (1141) leads to V q with modified 
quark mass and chemical potential. Note that the quark mass is also modified by 
Z x . Finally the effective potential is found to be, 



Fctt =-^"aux + V q (m g ; p,, T) , (15) 

^ =\v 2 + - «?) + y - N c log Z x , (16) 

rhq =(b a a + m )/Z x . (17) 



We find that the plaquette contributes to the effective potential in the modification 
of the wave function renormalization factor Z x , the quark mass m g , and the shift of 
the effective potential \x in addition to some auxiliary field terms, J-g, ux . 
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Since the 1/g 4 contributions have ambiguities in the present NLO treatment, 
we here compare the results in several ways of truncation. In the first treatment, 
abbreviated as NLO-A, ^" e g in Eq. (|15p is treated as it is, and we do not invoke any 
further approximations. In the second treatment (NLO-B), 0(1/ g) contributions 
in Z x and p are truncated, and simplified as Z x = 1 + /3 T tp T and p = p, — j3 T u T . 
In this treatment, we find that (p T and u> T couple to quarks separately through m q 
and p in V q , respectively. In the third prescription (NLO-C), we further truncate 
0(1/ g 4 ) terms in rhq and in logZ x . It is also possible to expand V q with respect to 
5fj, = fj, — fi, (NLO-D), 

V q (m q ; p, T) ~ V q (m q ; p, T) - f3 T uo T dV q /dp . (18) 

In each treatment NLO-A, B,C, and D, we evaluate J- e a under the stationary 
condition with respect to the auxiliary fields, <P = a,ip a ,tp T ,u} T , 

dFeft _ d-Faux dV q drhq dV q dp _ 
d$ ~ d$ drh q d$ dp d$ ~~ ' [ ' 

Note that V q depends on the auxiliary fields via the two dynamical variables fh q and 
p. Here we show the stationary conditions in NLO-A, as an example. Substituting a 
for ^ in Eq. (|19|) . we obtain the relation, a = —(l/Z x )(dV q /dm q ). By utilizing this 
result, the stationary condition for ip s leads to ip s = a 2 . We can solve the coupled 
equation for the stationary conditions of (p T and uj t as, 

2( A) dV q _ 

l Pr = — n ■ . 3 UJ T = -— = p q , 20 

where cpo = N c — Z x m q + /3 r u^. The second equation in Eq. ([20]) indicates the 
auxiliary field oo T is nothing but the quark number density p q . The equilibrium 
condition is determined self-consistently by minimizing T e s i n terms of a under 
the constraint lo t = p„(T, fj,;a,u T ). Stationary conditions in NLO-B, C and D are 
solved similarly. In NLO-D, uj t is explicitly obtained as a function of a, uj t = 
—dV q (m q (a);ii,T)/dp,, where the rhs does not contain uj t . In this meaning, the 
NLO-D gives a similar formulation to that in the previous work. 7 ) 

One of the most interesting features in the multi-order parameter treatments 
(NLO-A, B and C) is that it predicts the existence of partially chiral restored (PCR) 
matter at N c = 3 for relatively large j3 = Q/g 2 values. In Fig. [3l we show the phase 
diagram at (5 = 4.5 in the chiral limit. In NLO-A and B, the highest temperature 
of the first order phase boundary decreases, and the critical point deviates from the 
second order phase transition boundary at 6/g 2 — 4.5 and 3.0 in NLO-A and NLO-B, 
respectively. In NLO-C, the second order critical chemical potential p^ nd ^ at T = 
overtakes the first order one at 6/g 2 ~ 3.5. Between the first and second order phase 
boundaries, we find PCR matter. 

In Fig. U we show the comparison of p q and a in the present treatments NLO- 
A, B, C and D. The gradual increase of the quark density is a common feature 
of the multi-order parameter treatments. At low temperatures, we can investigate 
the appearance of the medium density matter more intuitively. The quark number 
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Fig. 3. Phase diagram in NLO-A, B, C and D with 7V C = 3,6/ g 2 = 4.5 Solid and dashed lines 
show the first and second order phase transition boundary, respectively. 

density p q = —dT^/dp is evaluated as, 



where x = exp[—(E q — p)/T\. When E q > jl is satisfied at small T, we obtain 
x — > and p q — > 0, while E q < jl leads to x — ► oo and — ► iV c . Medium density 
< p q < N c can appear only in the case where the energy and chemical potential 
balances, E q = jl, and x stays finite at T = 0. Since jl is a decreasing function of u T , 
we may have a medium density solution of Eq. (|20p in the region /2(<t, w t = iV c ) < 
E q (a,uj T ) < p. Specifically in NLO-B and C, E q = jl is found to be equivalent to 
the density condition p q = (p — E q )/(3 T , which can take the a medium value. In the 
large (3 region, this medium density matter can emerge in equilibrium as indicated 
in Fig. UJ Also in NLO-A, medium density matter appears in a similar mechanism. 
Thus the multi-order parameter treatment is essential to obtain the medium density 
matter at low T, and we observe chiral transitions twice as p increases. After the 
first one, the condition E q = jl is approximately satisfied. This state may correspond 
to the quarkyonic matter where we expect p ~ (constituent quark mass). 14 ) 

In summary, we have evaluated the effective potential in the next-to-leading or- 
der strong coupling lattice QCD (NLO SC-LQCD), where 1/g 2 effects are taken into 
account. We have discussed the chiral transition and the possibility of the quarky- 
onic transition. Here the order parameter of the latter is the quark number density, 
which can be naturally introduced via the NLO effects. The 0(1/ g 4 ) ambiguities 
have been examined by comparing the results in several truncation schemes, and we 
have found the following common properties as far as the quark number density is 
treated as the order parameter in addition to the chiral condensate. (I) The partially 
chiral restored (PCR) matter can appear in the large Q/g 2 region, (II) PCR sits next 
to the hadronic Nambu-Goldstone (NG) phase in the p direction, (III) the quark 
number density is high as 0{N C ) in PCR, (IV) after the "NG^PCR transition", 
the effective chemical potential is approximately the same as the quark excitation 
energy, (V) and the chiral transition to the Wigner phase follows after NG— >PCR 



2sinh[iV c /2/T] 



(21) 



N c X Nc + 2 cosh [N c ~p/T] t^o l + x N ° ' 
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Fig. 4. In the upper panel, solid and dashed curves show the first and second order phase transition 
boundaries, and dots show the critical end point. In the middle and lower panels, solid, dotted, 
dashed, and dot-dashed curves show the results in NLO-A, B, C and D, respectively, and dots 
and open squares show the first and second order transition points. 



transition. All these properties would be the essence of the quarkyonic matter and 
transition proposed in Ref. 1 ** In the previous work, the quark-driven Polyakov loop 
evaluated in SC-LQCD is shown to be small as 0{1/N C )^ and it would not grow 
much at low temperatures. This feature also agrees with the proposed property of the 
quarkyonic matter. In the present analyses, we have found the two sequential chiral 
transitions can occur along the /x direction at low T, and the first one involves the 
quarkyonic transition. Thus we can conclude that we have examined the quarkyonic 
picture with the clear connection to the finite coupling effects in the strong coupling 
expansion. The detailed analyses of the phase diagram evolution due to the finite 
coupling effects will be shown elsewhere. 16 ) 

There are many points to be improved in the present analysis. First, 0(l/g ) 
ambiguities is not small, and it is necessary to extend the analysis to the next-to- 
next-to-leading order (NNLO). Higher orders in 1/d expansion may be also necessary. 
The quarkyonic transition may be related to the so-called baryon mass puzzle; 9 ^' 
N c fj, c is calculated to be smaller than the baryon mass in the strong coupling limit. 
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This means that the chiral transition takes place before nuclear matter is formed. 
The chemical potential shift discussed in this paper may be the key to solve this 
problem. Competition with the color superconducting (CSC) phase and comparison 
with the MC results at finite baryon density 10 )' 15 - 1 are other interesting subjects to 
be investigated. 
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